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Abstract 


The most distressing types of heart malfunction occur because of an abnormal rhythm of 
the heart. Cardiac arrythmias can be caused by abnormal rhythmicity of the pacemaker, 
electrolyte disturbances, blockage of the transmission of the electric impulse through the 
heart, and other abnormalities. There is strong evidence that space flight is associated 
with decreased cardiac electrical stability that may pose a life threatening risk to 
astronauts. For example, during the Skylab missions, a crewmember had a five beat run 
of ventricular tachycardia during lower body negative pressure. Also, analysis of nine 
24-hour Holter monitor recordings obtained during long term spaceflight on Mir revealed 
one 14-beat run of ventricular tachycardia. A Mir cosmonaut was replaced in 1986 
because of cardiac dysrhythmias. Most recently, in July of 1997, a Mir commander was 
unable to participate in the Spektr module repair due to complaints of an irregular heart 
rhythm. Despite these examples, possible mechanisms of arrhythmias and 
countermeasure strategies have barely been addressed. 

The Poincare method has been proposed as a technique that might potentially predict life- 
threatening arrhythmias before they occur. According to this method, each RR interval 
obtained from an EKG recording is plotted sequentially vs. the previous RR interval. 
Several studies using the method have demonstrated a strong correlation between the 
shape of the Poincare plot and ventricular arrhythmia. Our purpose was to develop an 
automated software program that detects the “R” peaks from an EKG recording while 
simultaneously displaying the Poincare plot and other related parameters. 
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INTRODUCTION 


The pumping action of a heart is regulated by a conduction system that spreads electrical 
impulses from the sinoatrial node to the rest of the heart. The signal generated in the 
sinoatrial node first disperses through the atrial tissues causing the atria to contract. The 
signal then moves through the Purkinje fibers into the ventricular septum. Finally, the 
signal is conducted to the outer ventricular walls causing the ventricles to contract. In a 
normal heart, this process is very organized. However, in a diseased heart, the 
conduction system can be damaged. The degree of irregularity in the conduction of the 
impulse depends greatly on the extent and location of the disease in the heart. 

Electrical potentials that originate in the heart spread out into neighboring tissues. 
Electrodes that are placed on the chest opposite to the heart can record these potentials. 
Such a recording is known as an electrocardiogram (EKG). A sample EKG is shown in 
Figure 1 . There are three distinct peaks in an EKG signal: the P wave, the QRS complex, 
and the T wave. Electrical potentials, generated when the atria depolarize at the onset of 
atrial contraction, cause the P wave. The QRS complex is caused by potentials generated 
by ventricular depolarization before contraction. When the ventricle recovers from 
depolarization, the T wave is generated. 



Figure 1 . A Sample EKG Showing the Characteristic P, Q, R, S and T Waves. 


The P-Q or P-R interval is the amount of time it takes from the onset of atrial contraction 
to the onset of ventricular contraction. The Q-T interval is the time of ventricular 
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contraction. The time between successive QRS complexes is the time it takes for a 
complete heart beat. This is also known as the RR interval. 

By analyzing the peaks and intervals in an EKG, clinicians can determine whether the 
heart is functioning properly or if it is diseased. Such analysis also provides a good way 
of determining the type of cardiac disease that may be present. Recently, there has been 
an increasing interest in one method that relies on the RR interval. This method, called 
the Poincare plot (2,5,12), is simply a plot of the current RR interval (RRn+i) vs. the 
previous RR interval (RRn). If A, B, and C are consecutive RR intervals, their Poincare 
plot consists of two points with the coordinates (A,B) and (B,C). This is shown in Figure 
2. When EKG recordings are of longer duration (i.e., minutes to hours), the Poincare 
method provides a visual way of analyzing overall heart rate variability (HRV) (5,8). 
When the RR intervals are very similar, i.e., low standard deviation, the points in the 
Poincare plot form a tight pattern. When the RR intervals vary a lot from beat to beat, 
i.e., high standard deviation, the points on the Poincare plot are spread out. Since several 
studies suggest that both heart failure and impending arrhythmias can be associated with 
a decreased HRV (3,4,7), the Poincare plot provides a useful diagnostic tool. 


Poincare Plot 




Figure 2. Poincare Plot of Three Consecutive RR Intervals 


Healthy subjects usually have a so called “comet” shaped (13) 24-hour Poincare plot 
(Figure 3). The range of RR intervals is usually between 500 to 1000 ms. On the other 
hand, 24-hour Poincare plots obtained from patients with heart failure may exhibit three 
different abnormal patterns. The most common one is the “cigar” shape. This shape is 
attributed to little or no increase in R-R variability at lower heart rates. The second 
Poincare plot pattern from unhealthy subjects is the “fan” shape. Both of the cigar and 
fan patterns have a much lower RR interval range. The third pattern is a “complex” 
pattern that is a combination of the cigar and fan patterns that are spread out over the 
entire plot. 


20-4 





Comet 


Cigar 




Fan 

Complex 

A 


V ’C 



* i 





Figure 3. Different Patterns of a 24-hour Poincare Plot. 


METHODS 

The first step in generating a Poincare plot is to acquire the locations of the R-peaks in 
the signal. To accomplish this goal, an R-wave detector was developed that functions as 
follows: First, the Fast Fourier Transform (FFT) of the EKG signal is determined. Next, 
the lower frequencies of the Fourier spectrum are filtered out. Figure 4 shows the Fourier 
spectrum of an EKG signal before and after filtering. By filtering the spectrum, the 
slower components of the EKG signal are removed. Finally, the Inverse Fast Fourier 
Transform (IFFT) is computed from the filtered spectrum, and a modified version of the 
EKG signal is obtained. The only significant peaks in this modified signal are the R 
waves (Figure 5). 

The locations of the R peaks were determined using a custom made peak detector. This 
detector first divides the signal up into different windows. For each window it then 
computes the maximum. Finally, the detector scans through the maxima to determine 
which windows have a maximum that has a higher value than the two maxima in the 
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surrounding windows. The location of that maximum in the EKG signal is then 
determined and saved as an R peak. 


A 



Figure 4. (A) FFT of EKG Signal Before Filtering. (B) FFT After 

Filtering. 



Figure 5. The EKG Signal Before (solid) and After (dashed) Filtering. 
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From the locations of the R peaks, the RR intervals can then be calculated. Most of the 
important information regarding HRV can be obtained from the RR intervals. The 
following is a short list of some of these variables and their descriptions (9). 


RRmcan 

Average of RR intervals 

RRstd 

Standard deviation of RR intervals 

ARR 

Difference between elements of RRn+i and RRn 

ARRstd 

Standard deviation of ARR 

RR„std 

Standard deviation of RR n 

RRn+istd 

Standard deviation of RRn+i 

nRR n std 

RRnStd / RRstd 

nRR n +istd 

RRn+istd / RRstd 

NN50 

Sum of elements of ARR that are less than 50 ms 

pNN50 

NN50 divided by length of ARR 


By plotting RR«+i vs. RRn, a Poincare plot is generated. Once the Poincare plot is 
acquired, an ellipse (3) is fitted to the data (Figure 6). To produce such an ellipse, the 
confidence interval method was used. Knowing the means (pi,p2), standard deviations 
(cri,<T 2 ) and the correlation coefficient (p) of RR n and RRn+i, one can define an ellipse that 
encloses a desired percentage of the data. The equation for such an ellipse is the 
following: 


1 

2(1 -P 2 ) 




2p(x-/u x )(y-ju 2 ) | {y-Mi) 
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Eq. 1 


where n is the desired percentage of data to be enclosed within the ellipse. In Equation 1, 
one can assume that RRn and RRn+1 have the same mean (pi = p. 2 = p) and standard 
deviations (a i = 02 = a), which simplifies the equation into the following: 
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Poincare Plot 
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Figure 6. A Poincare Plot with its Ellipse and Axes. 


The only variables in this equation are x and y. Therefore, by substituting for either of 
these variables, the other variable can be determined. Using x and y the ellipse can be 
superimposed on the Poincare plot. The major and minor axes of the ellipse can then be 
determined. These axes are directly proportional to RR„std and RRn+istd, respectively. 

The histograms of RR and ARR can also be determined (Figure 7). These histograms are 
then analyzed to calculate the HRV triangular index and the triangular index of RR 
(TIRR) (7). The HRV triangular index is defined as the number of RR intervals divided 
by the maximum of the RR histogram. TIRR is the base (B) of the triangule fit to the 
histogram. To find B, the area under the histogram is first computed. The area of a 
triangle is the product of B and the height (H), divided by 2. If H is equal to the 
maximum of the histogram, then B can be calculated from the area of the triangle. To 
plot the triangular fit, the ratio of the area of the histogram before the maximum to the 
area of the histogram after the maximum is first calculated. The base is then positioned 
so that the ratio of its length before the maximum to its length after the maximum is equal 
to the ratio of the areas. 
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Figure 7. RR Histogram Showing the Triangular Fit. 


RESULTS 


Software Package 

To carry out these methods and calculations in an automatic an easy manner, a graphical 
user interface package was developed using Matlab. This package allows the user to 
open several formats of EKG signals. The electrocardiogram is immediately displayed 
when the data file is successfully opened. The user can then analyze the data using 
several options in the package. The first option is detecting the R peaks in the signal. 
When this method is called from the menu, a window like that shown in Figure 8, is 
opened. In this window, there is a plot of the EKG signal with the R-wave peaks marked 
by asterisks, a slider to scroll through the signal, an edit box to change the size of the 
window used to detect the R peak, and a [Recalculate] button that recalculates the peaks 
according to the window size that is defined in the edit box. 

The second option in this package is the Poincare analysis. When this menu item is 
chosen, the Poincare plot is displayed along with the ellipse that fits a default percentage 
of data. This plot also shows the major and minor axes. The third option is to plot the 
histograms of RR and ARR. 
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Figure 8. The R Detector Window 


It is also possible to generate a 3D Poincare plot (10,1 1). In such a plot (Figure 9), each 
point’s coordinates are (RR n +2, RR n +i, RR n ). Once the analysis is finished, the user can 
print out a list of all the variables associated with the methods described in this paper. 



Figure 9. 3D Poincare Plot Example 
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Data Analysis 

Two EKG signals obtained from a healthy and non-healthy subject were analyzed to 
obtain the Poincare plots shown in Figure 10. There is a remarkable difference between 
the healthy subject’s Poincare plot and that of the unhealthy subject. This shows that the 
methods integrated into this software provide reliable data. 


♦ A B 




Figure 10. Poincare Plots From a (A) Healthy Subject and a Subject with 
Heart Disease 


CONCLUSION 

In this project, the Poincare method of analyzing HRV was automated. The software 
package that was developed allows the user to open an EKG file and process the data to 
obtain a Poincare plot, histograms of RR intervals, and variables that can help physicians 
determine whether the patient’s heart is functioning normally or abnormally. 

The Poincare package developed in this project can only be used in a Matlab compiler. 
Therefore, for further research, the methods involved in this package should be integrated 
into a standalone software package. The R wave detector is also a very rudimentary one. 
A more sophisticated R wave detector where no user input is needed would automate the 
software even further. Also, this package works offline only. It would be a greater asset 
for physicians to be able to monitor the Poincare plot in real time. 

Several other analysis tools can be added to this package to enhance its diagnostic use. 
The approximate entropy (ApEn) method (1,6) for assessing complexity of the EKG 
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signal provides a sensitive marker for changes in HRV. The power spectral density 
(PSD) of RR intervals can also be added to this package. PSD (2,9) is used to study the 
sympathetic and parasympathetic neural control to the heart. 
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